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Abstract 

We analyse the anomalous properties of specific electronic states 
in the Kronig-Penney model with weak compositional and structural 
disorder. Using the Hamiltonian map approach, we show that the 
localisation length of the electronic states exhibits a resonant effect 
close to the band centre and anomalous scaling at the band edges. 
These anomalies are akin to the corresponding ones found in the An- 
derson model with diagonal disorder. We also discuss how specific 
cross-correlations between compositional and structural disorder can 
generate an anomalously localised state near the middle of the energy 
band. The tails of this state decay with the same stretched-exponential 
law which characterises the band-centre state in the Anderson model 
with purely off-diagonal disorder. 

Pacs: 73.20. Jc, 73.20.Fz, 71.23.An 

1 Introduction 

The interest in one-dimensional random models has been constantly increas- 
ing since the discovery that, contrary to previous beliefs, systems of this class 
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can undergo a sort of localisation-delocalisation transition when disorder ex- 
hibits specific long-range correlations [JJ. Nowadays, one-dimensional models 
with correlated disorder are used in many fields of physics, including semicon- 
ductor superlattices [2], bilayered media [3], Bose-Einstein condensates [I], 
transmission in waveguides |5j and structures with corrugated surfaces [6], 
and DNA modelling [7J. 

In many of these problems, the system under study is represented in terms 
of an aperiodic Kronig-Penney model. The original Kronig-Penney model 
was introduced in the early 1930s to analyse the band structure of crystalline 
materials [8] and found new applications in the 1980s when it was used to 
study electronic states in semiconductor superlattices [SJ- More recently, the 
interest in random one- dimensional systems has spurred investigations of the 
aperiodic variants of the Kronig-Penney model, which have been analysed 
both theoretically and experimentally [2J (5J [TO] . 

Because of the simplicity and versatility of the aperiodic Kronig-Penney 
model, it is important to obtain as much analytical information as possible on 
the structure of its electronic states. A detailed study of the Kronig-Penney 
model with weak compositional and structural disorder was done in [TTj . 
where an analytical expression for the localisation length of the electronic 
states was derived. The formula works very well for most values of the 
energy, but fails in a neighbourhood of the band centre and of the band edges, 
where localisation anomalies appear. The first goal of the present work is to 
understand these anomalies: using the Hamiltonian map approach [12J, we 
show that they have the same form of the corresponding anomalies which 
exist in the standard Anderson model with diagonal disorder at the band 
centre [TBI [Ti] and at the band edge [15]. These anomalies have recently 
attracted attention because of the related violations of the single-parameter 
scaling hypothesis [16J at the edge [T7] and at the centre [TS] of the energy 
band. It is therefore of some interest to see whether similar anomalies exist 
in the random Kronig-Penney model. 

After explaining the nature of the anomalies of the localisation length 
at the band centre and band edges, we focus our attention on a particular 
variant of the random Kronig-Penney model, characterised by specific cross- 
correlations between the compositional disorder and the structural one. We 
show that such a Kronig-Penney model corresponds to a Anderson model 
with energy-dependent diagonal and off-diagonal disorder. For a specific 
value of the energy, the diagonal disorder vanishes and the Kronig-Penney 
model becomes equivalent to the Anderson model with purely off-diagonal 
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disorder. As is known [HI], at the band centre the Anderson model with 
purely diagonal disorder has anomalously localised states, whose existence 
is another violation of the single-parameter scaling hypothesis (20]. Because 
of the identity of the mathematical equations, similar states with stretched- 
exponential tails also exist in the considered variant of the Kronig-Penney 
model. 

This paper is organised as follows. In Sec. [2] we define the model under 
study and we summarise the main results concerning the localisation of the 
electronic states in the standard case. In Sec. [3] we discuss the anomalies of 
the localisation length emerging for values of the energy close to the band 
centre and in the neighbourhood of the band edge. In Sec. H]we analyse the 
anomalously localised states created by particular cross-correlations of the 
disorder. The conclusions are drawn in Sec. |5j 

2 Localisation of the electronic states: the 
standard case 

In this section we define the model under study and we briefly derive the 
expression for the localisation length valid in the general case. 

2.1 Definition of the model 

We consider a Kronig-Penney model with weak compositional and structural 
disorder. The model is defined by the Schrodinger equation 

— ip"{x) + (U + u n )5{x — an — a n )ip(x) = q 2 ip{x) (1) 

n=— oo 

which describes the motion of an electron of energy q 2 in a potential formed 
by a succession of aperiodically positioned delta-barriers of random strengths. 
To simplify the form of the equations, here and in what follows we use energy 
units such that h 2 /2m = 1. We introduce compositional disorder in model ([1]) 
by assuming that the barrier strengths display random fluctuations u n around 
their mean value U. The structural disorder, on the other hand, is present 
because the positions of the barriers are displaced by a random amount a n 
with respect to the lattice sites na. 
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For weak disorder, it is enough to specify the statistical properties of 
the model in terms of the first two moments of the strength fluctuations u n 
and of the relative displacements A n = a n +i — a-n (which are more relevant 
than the absolute displacements a n themselves for the description of the 
electronic states). We assume that both variables have zero average, (u n ) = 
and (A n ) = 0, and that their variances satisfy the conditions of weak 
compositional disorder 

«) « u\ 

and of weak structural disorder 

( A n)<? 2 « 1, {^l)U 2 « 1. (2) 

We also consider the normalised binary correlators 

m (u„u n+ i) 

xM = ~W~ 

(u n A n+ i) 



Xa(0 



(u n A n ) 



as given functions. We do not attribute any specific form to the correl- 
ators (E]); we only assume that, because of the spatial homogeneity in the 
mean of the model, they are (decreasing) functions of the difference / of the 
site indices. For the sake of simplicity, we restrict our attention to the case 
in which the self-correlators an d are even functions of I, but we 

do not make the same assumption about the cross-correlator X3(0 (i n other 
words, we suppose that the model is only partially isotropic). 



2.2 The Hamiltonian map approach 

The Schrodinger equation (pQ) has the same form of the dynamical equation 
of a stochastic oscillator with frequency q perturbed by a random succession 
of delta-shaped impulses ("kicks"). The Hamiltonian of such an oscillator is 



CO 

E 



(U + u n ) 5 (t — an — a n ) 



x 



(4) 
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The mathematical identity of the dynamical equation of the kicked oscil- 
lator PJ with Eq. (JTJ) allows one to analyse the properties of the elec- 
tronic states of the Kronig-Penney model ([1]) in terms of the trajectories 
of the kicked oscillator (j3J): this is the key idea of the Hamiltonian map ap- 
proach [T2l [2T] . After integrating the Hamiltonian equations of the stochastic 
oscillator over the interval between two successive kicks, one obtains the 
Hamiltonian map 



x n+l 
Pn+1 



Pn 



(5) 



with the transfer matrix 



cos [q (a + A n )] + (U + u r 
-gsin [q (a + A n )] + (U 



- sin [q (a + A n )] - sin [q (a + A n )] 
q q 

u n ) cos [q (a + A n )] cos [q (a + A n )] 

(6) 

To analyse the trajectories of the Hamiltonian map flSJ), we follow the ap- 
proach proposed in [TU] (see [TTJ for details) and we perform the canonical 
transformation (x n ,p n ) — > (X n ,P n ) of the form 



x n 

Pn 



( 



\ 



a cos ■ 



qa 



-qa sin 



qa 



1 qa \ 

— sin — » 

qa 2 
1 qa 

— cos — 
a 2 



(7) 



The parameter a in Eq. (jjTJ) is defined by the identity 

U 



^ sin (qa) 



a: 



2q 



[cos (qa) — 1] 



q 2 U 

sin (qa) — — [cos (qa) + 1] 
Zq 



The utility of the transformation (|7J) lies in the fact that the Hamiltonian 
map (jSJ), when expressed in terms of the new variables (X n ,P n ), is reduced 
to a simple phase-space rotation in the absence of disorder. The rotation 
angle ka of the unperturbed map is defined by the identity 



cos (ka) = cos (qa) + — sin (qa) . 

2q 



(8) 



Eq. OH]) determines the band structure of the Kronig-Penney model (pQ) and 
reveals that the parameter k is the Bloch wavenumber. 
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To simplify the form of the Hamiltonian map fl5]) further, we switch from 
the Cartesian coordinates (X n , P n ) to action-angle variables ( J n , 6 n ) via the 
standard relations 

X n = a/2 J n sin 6 n 
P n = \/2J^cos9 n ' 

Within the second-order approximation, the Hamiltonian map fl5]) can be 
written in terms of the action-angle variables in the form [TT] 



7 _ f)2 7 

"ri+l -^n^Ti 

1 1 

6 n +i = On + ka - - [1 - cos (26 n + fca)] u n + - [C - cos (26 n + 2ka)\ A r 
1 

+ - [2 sin (26 n + ka) - sin (40 n + 2fca)] u 2 

+ - [2C sin (20 n + 2ka) - sin (4fl n + 4ka)] A 2 n 

+ — [sin (ka) — 2 sin (2# n + 2ka) + sin (4# n + 3ka)] u n A n 

with the ratio of the action variables being equal to 



(9) 



Dl = 1 + sin (26 n + ka) u n - sin (26 n + 2A;a) A n 

+ - [1 - cos (20 n + ka)} u 2 n + i [1 - C cos (20 n + 2ka)] A 2 
+ [cos (26 n + 2fca) — cos (ka)} u n A n . 

In Eqs. 09]) and ffTOl) we have introduced the rescaled disorder variables 

sin (qa) 



(10) 



g sin (ka 



-u 



,., and A n 



U . 
sm (reaj 



and the parameter 



c 



q sin (fca) 

c7 



ga 2 + 



qa' 



(11) 



(12) 



2.3 The localisation length 

The inverse localisation length for the Kronig- Penney model ([1]) is defined as 



r 1 



loc 



lim , 

7v->oo IS a 



1 N 



71+1 



71 = 1 
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In the dynamical picture, Z lo * is equivalent to the Lyapunov exponent of the 
Hamiltonian map 



A 



1 N 
lim > lo 

iv^oo Na ^ 



n=l 



■'n+l 



^ E lo s ^ £ lo s (*») > 

(13) 



Af->oo A^a ^— ' ° v ""' iV-+oo Aa 

n=l n=l 



with D n being defined by Eq. ( flQl) and 

>/C + 1 [1 + cos (ga)] sin 6 n+1 



\JC, — 1 sin (ga) cos 



a/C + 1 [1 + cos (ga)] sin 6 n + a/C — 1 sin (ga) cos 9 n 



(14) 



As discussed in Sec. [31 in the weak-disorder case the invariant distribution for 
the angular variable is either flat (the normal case) or presents a moderate 
modulation (as happens at the band centre and for the other cases in which 
the Bloch wavenumber is a rational multiple of n/a). Only at the band 
edge the invariant distribution is strongly modulated even for weak disorder. 
Therefore as a function of n the ratio ( fl4l) can be expected to oscillate around 
a unitary value, unlike the ratio D n which is larger than one on average 
because of the exponential increase of the action variable. As a consequence, 
away from the band edge, one can drop the second term in the right-hand 
side of Eq. (113"]) and compute the inverse localisation length as 



A 



1 N 
7V->-oo A a z — ' 



(15) 



n=l 



Within the second-order approximation one can expand the logarithm in 
Eq. (TT5T) and obtain 



A 



+ 
+ 



sin (26 n + ka) u n — sin (26 n + 2ka) A r 



[1-2 cos (20 n + ka) + cos (40 n + 2ka)\ u 2 n 
[1 - 2Q cos (26 n + 2ka) + cos (46 n + 4ka)] A r 2 



(16) 



— - [cos (ka) — 2 cos (26 n + 2ka) + cos (A6 n + 3ka)] w n A n 



In the general case, the average over the angle variable in the right-hand side 
of Eq. f fl6|) can be computed using a flat distribution p(6) = l/(2n). In fact, 
as the Hamiltonian map (jH]) shows, the angle variable has a fast dynamics 



7 



compared to the action variable and quickly assumes a uniform distribution 
in the interval [0, 2tt]. In this way, following the steps of Ref. [UJ, one obtains 
that the inverse localisation length (j!5p has the form 



A 



1 



(vftWtika) + (A 2 n )W 2 (ka) - 2(u n A n )W 3 (ka) 



(17) 



8a L 



with 



Wxika) 



Xi{l) cos (2kal) , 



l=— oo 
+oo 



W 2 (ka) 




l=— oo 
+oo 



W 3 (ka) 



X3(l) cos [ka (21 + 1)]. 



Eq. (iT7j) constitutes the standard expression of the inverse localisation length 
of the electronic states in the Kronig- Penney model (0Q). Note that expres- 
sion (117j) is slightly more general than the corresponding formula given in [TT] 
because we have dropped the assumption that the cross-correlator Xs(k) is 
an even function of k. The difference between the two expressions shows up 
in the last term of Eq. f|T71) and in the definition of the W 3 (ka) power spec- 
trum. When Xz(k) — X3(~^), Eq. ffTTj) reduces to the inverse localisation 
length derived in 

The inverse localisation length (|17p works well for almost all values of 
the energy, but fails close to the band centre and at the band edges, where 
anomalous effects appear. The nature of these anomalies is discussed in the 
next section. 

3 The band-centre and band-edge anomalies 
of the localisation length 

As mentioned in the previous section, the derivation of formula (TTT1) rests on 
the crucial assumption that the angle variable of the map (Q has a uniform 
invariant distribution. This assumption is generally justified on the grounds 
that the angle variable evolves much faster than the action variable and 
quickly sweeps the whole interval [0, 2ir], One should notice, however, that 
this argument cannot be applied when the Bloch wavevector is a rational 
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multiple of ir/a, because in this case the noiseless angular map has periodic 
orbits, whose effect persists in the form of a modulation of the invariant 
measure when a weak noise is switched on . For weak disorder, and within 
the limits of the second-order approximation, this implies that the general 
formula ( ITT)) cannot be applied when the Bloch vector takes values ka ~ 
or ka ~ ±7r, i.e., at the edges of the energy band, and for k = ±7r/2a, i.e., 
for the energy [q(ir /2a)] 2 which lies close to the band centre. In principle, 
one should expect anomalies for all Bloch vectors of the form ka = n/n 
with |n| > 2. In practice, however, in these cases the invariant distribution 
of the angular variable is modified by a perturbative term proportional to 
cos(2n#) or sin(2n#) and this does not affect the outcome of the average in 
Eq. fllSp . which contains only second- and fourth-order harmonics of 9 (the 
same conclusion applies to the Anderson model, as discussed in Ref. |14j). 

In conclusion, within the second-order approximation anomalies are found 
only at the band edge {ka ~ or ka ~ ±7r) or close to the middle of the 
band [ka = ir/ '2). For these special values of k the assumption of a flat 
invariant distribution must be abandoned and the specific form of p(6) has 
to be determined before one can compute the localisation length of the elec- 
tronic states. The non-uniform distribution of the angular variable produces 
deviations from the standard formula (flT)) of the inverse localisation length; 
this section is devoted to the analysis of these anomalies. 

For the sake of simplicity, in this section we will restrict our attention 
to the case of uncorrelated disorder. In other words, we will consider two 
successions of random variables {u n } and {A n } such that 



(u n u k ) = (ul)S nk , (A n A fc ) = (Al)S nk , (u n A k ) = 0. 



In this case the standard expression (H 
reduces to the simple form 



of the inverse localisation length 



A 



1 

8a" 



ul) + (K) 



8a sin 2 (fca) 



sin 2 (ga) 



u 2 n ) + U 2 (Al) 



18) 



3.1 The anomaly near the middle of the energy band 

We will first consider the anomaly for ka = ±7r/2. The corresponding energy 
lies close to middle of the band and therefore we will often speak, somewhat 
loosely, of band-centre anomaly. By assigning the value k = ir/2a to the 
Bloch vector and by taking into account that the noise-angle correlators 
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vanish for uncorrelated disorder, one can reduce the general expression ( 1161) 
to the simpler form 



1 



A = — {[l + 2(sin(20 n ))-(cos(40 n )>] 



ul) 



<« 1 (19) 
+ [1 + 2C(cos(20 n )> + (cos(40 n ))] (A*)] . 

To determine the form of the invariant distribution p(8), we follow the 
method introduced in [12] and we consider the map for the angle variable in 
Eq. (Q which, for uncorrelated disorder and k = it /2a, simplifies to 

n+1 = 9 n + £ - - [1 + sin (20„)] u n + - [C + cos (20„)] (A n ) 

1 1 

+ - [2 cos (26 n ) + sin (40 n )] (u 2 n ) - - [2( sin (20 n ) + sin (40 n )] ( A 2 n ) . 

o o 

To get rid of the constant drift term, one can consider the fourth iterate of 
this map, i.e., 

n+ , = 9 n +±sm(49 n )((u 2 n )-(A 2 n ) 

- -[1 + sin (29 n )} (u n + u n+2 ) - - [1 - sin (29 n )} (u n+1 + u n+3 ) 

Going to the continuum limit, one can replace this map with the Ito stochastic 
differential equation 



{ -f-[l + sin (29)} dWi - J K -f- [1 - sin (20)] d^ 2 (20 ) 



+ \/ ^ [C + cos (20)] d^ 3 + y ^ [C - cos (20)] dw 4 

where Wi(t), . . . , W±{t) represent four independent Wiener processes. Given 
an initial condition 9(t ) = O , the Ito equation (1201) defines a stochastic 
process 9(t), whose conditional probability p(9,t\9o,t ) = p can be obtained 
by solving the associated Fokker-Planck equation 



2 



i d 

+ 4902 



3(u 2 n ) + (2C 2 + 1) <A*>) + ((A 2 ) - <«*>) cos (40)] p) 



(21) 
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with the initial condition p(9,t \9 ,t ) = 8(9 — 9 ) [22]. The stationary solu- 
tion of Eq. ( 121 p which satisfies the conditions of normalisation and periodicity 
is pj] 

y/A+\B\ 1 



(22) 



4K(C) ^A-Bcos {46)' 

In Eq. ( 1221) K(C) is the complete elliptic integral of the first kind and we 
have introduced the constants 



A = 3(^) + (2C 2 + 1) (A 2 ), 



^ = «> - (A„) 



and 



C 



2\B 




A+\ 


B\ 



\ 



2 


«) ~ (A 2 ) 






(ul) - {AD 


+ 3{uD + (2C 2 + 1) <A»> 



(23) 



Eq. ( 1221) shows that, as expected, when the Bloch vector takes the value k = 
7i /2a, the invariant measure has period tt/2. The numerical computations 
agree well with formula ( 1221) as can be seen in Fig. CD The data represented 
in Fig. [I] were obtained for mean field U — 8 and disorder strengths ■>/ {uD = 
a/ (A^) = 0.02. Here, and in the rest of the paper, we present numerical data 
which were obtained for energy values within the first band. 

The knowledge of the invariant distribution (T22]) makes possible to com- 
pute the Lyapunov exponent ( JT9|) . The averages of the functions of argument 
29 vanish, but the the average of cos(46 l ) does not and gives rise to the an- 
omaly of the localisation length. After some algebra, one obtains that the 
inverse localisation length for ka = tt/2 is 



A = — 



ul) - {K) 

ul) + C{K) 



3(^)+(2C 2 + l) {AD 



E(C) 
K(C) 



(24) 



where E(C) is the complete elliptic integral of the second kind and the ar- 
gument C is defined by Eq. (T23l) . The numerical computations confirm the 
existence of an anomaly for ka = tt/2 as can be seen from the data rep- 
resented in Fig. [2] which show a small but clear deviation from the value 
of the localisation length predicted by the standard formula ( fTHl) . The nu- 
merically computed inverse localisation length for ka = tt/2, on the other 
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0.19 
0.18 
0.17 
9 0.16 
0.15 
0.14 
0.13 

-K -7C/2 7C/2 K 

e 

Figure 1: (Colour on line) Invariant distribution p versus 9. The solid line 
corresponds to the prediction of Eq. (|22l) . while points represent numerical 
results. 



hand, matches well the theoretical value f )2^|) . We observe that the numer- 
ical data represented in Fig. [2] were obtained for a specific realisation of the 
disorder; when different disorder realisations are considered, the discrepancy 
between the numerical value of the Lyapunov exponent for ka = 7r/2 and the 
predicted result ( 124)) fluctuates slightly around zero, always assuming small 
values as in Fig. [2J 

In conclusion, in the Kronig-Penney model ([Q) a resonance effect occurs 
for k = 7r/2a and produces an anomaly of the localisation length. The effect 
has the same nature of the band-centre anomaly found in the Anderson model 
with diagonal disorder [T3(, [12] . 

We remark that, although our analytical results are restricted to the 
case of uncorrelated disorder, we found numerical evidence that disorder 
correlations can enhance the anomaly of the localisation length near the 
middle of the energy band, in agreement with the theoretical conclusions 
of Ref. |23j. As an example, we can consider the case of structural and 
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Theoretical prediction 
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(3) 
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k 0.0034 
\A 0.0033 
\ 0.0032 








■ 


0.007 


_ \ 0.0031 
\ (10030 












0.006 


\ 27T/5 


n/2 


37T/5 




0.005 






ka 








0.004 














0.003 














0.002 















7T/5 3^/10 2n/5 n/2 3n/5 7tt/10 4tt/5 

ka 



Figure 2: (Colour on line) Inverse localisation length A versus ka (uncorrel- 
ated disorder). The dashed line (1) represents numerically obtained values; 
the solid line (2) corresponds to formula flT8|) ; the symbol (3) represents the 
anomalous value fl24l) . The inset shows a close-up of the anomaly. 

compositional disorder with self-correlations of the form 

^ 3irl \5 / 

and no cross-correlations, XsiO = 0- The long-range correlations of the 
form (T25]) create mobility edges at k = n/ba and k = 4ir/5a (see [11] for 
details). Fig. [3] represents the numerical data obtained for this kind of dis- 
order with U = 8 and yjv^j = a/ (A^) = 0.02. The data clearly show an 
enhanced anomaly at k = 7r/2a with respect to the case of totally uncorrel- 
ated disorder. Adding cross-correlations does not introduce any significant 
modification to the picture. 

We would like to stress that, although the previous example shows how 
correlations of the disorder can enhance the anomaly for k = it /2a, not all 
correlations produce the same effect. This can be appreciated in the case 
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0.03 



0.025 
0.02 



A 0.015 
0.01 
0.005 


71/10 K/5 37T/10 27T/5 K/2 37T/5 77T/10 4k/5 9k/\0 k 

ka 

Figure 3: (Colour on line) Inverse localisation length A versus ka (self- 
correlated disorder). The dashed line (1) represents numerically computed 
data, while the solid line (2) corresponds to the values predicted by Eq. (fT7|) . 
The inset shows the anomaly in detail. 



analysed in Sec. HI in which the resonance effect for k = ir/2a is shadowed 
by a different kind of anomaly generated by specific cross-correlations for a 
value of k close, but not identical, to tt /2a. 

3.2 The band-edge anomaly 

We now turn our attention to the anomaly for ka = e — >■ + , i.e., in the 
neighbourhood of the band edge. Because of the similarity with the band- 
edge anomaly in the Anderson model, we can apply the method used in [12] 
to the present case. We first derive the form taken by the map (Q at the 
band edge. For ka = e — > + the rescaled random variables ( ITTjl can be 
approximated as 

sin(g a)u n ~ A n 
u n = h • • • and A n = U h . . . , 

q Q e e 
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where we have introduced the symbol go — Q(k = 0). In the same limit, the 
parameter (lT2"j) reduces to 



4g 2 



( = \I 1 + Jp s ^(ka) = l + T %e ' + ... 



U 2 



(26) 



Taking into account these approximations, one can write the Hamiltonian 
map in the form 



Jn+l 
9 n +l 



D 2 J 

9 n + e - — sin 2 ( 

e 



with 



D, 



1 sin 

e 



+ sm d 
e 2 



a 



COS I 



+ . . . 



cos [9 n ) + — sim 
le 2 



+ 



(27) 



(28) 



In Eqs. (1271) and (1281 . the symbol represents the linear combination of 
structural and compositional disorder 

t tt\ sin(q a) 



with zero average, (£ n ) = 0, and variance 



sin 2 (g a) ., 



Going to the continuum limit, one can replace the angular map in Eq. (J2 
with the stochastic Ito equation 



d9 



a 



e + — sin 3 (0) cos (0) 



a' 



dt + sin 2 (9) dW 



whose associated Fokker-Planck equation 



dp 
di 



d_ 

89 



a 



e + — sin d (0) cos(0) 



1 d 2 
P > + 2W 



a 



— sin (9)p 



(29) 



(30) 



gives the conditional probability p(9,t\9 ,t ) = p for the stochastic process 
9(t) [22J. By introducing the rescaled time 



a 

t = —t, 

e 2 
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one can cast the Fokker-Planck equation (1301) in the form 



dp 







1 d 2 



89 



2d9 2 



a 



— sin (6)p 



(31) 



which contains the noise intensity a 2 and the distance from the band edge e 
combined in the single scaling parameter 



x 



a- 



The invariant distribution p{9) is the stationary solution of the Fokker- 
Planck fl3T|) which is normalisable and satisfies the periodicity condition 
p{9 + 7r) = p{0). The solution possessing these features is [T5l [12] 



-/(*) 



N{x)sm 2 (6)J e sin 2 (0) 



=>/(*) 



(32) 



with 



and 



f(6) = 2k 



cot 3 (#) + cot(0) 



iV( 




2tt 



exp 



, i x 

-2k V x 

12 



dx. 



The integral representation (132]) defines the invariant measure in the inter- 
val [0,7r]; p(0) can be extended outside of this interval via the periodicity 
condition p{6 + ir) = p(9). 

To obtain a qualitative understanding of the behaviour of the invariant 
distribution (|32|) . it is useful to consider its values at the edges and at the 
centre of the [0, 7r] interval. For 9 — > + and 9 — > n~ one has 



5/3 



if k ->• 



— if x — )■ oo 

2tt 



while for # — > 7r/2 the invariant distribution behaves as 



P0) 



x 



1/3 if x -»■ 



-1 if 

2tt 



x — > oo 
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These equations show that, when the energy moves closer to the band edge 
on the scale defined by the disorder strength, i.e., for k — > 0, the invariant 
distribution develops two pronounced maxima for 9 ~ and 9 ~ n. This 
conclusion is supported by direct numerical computation of the invariant 
distribution, as shown by Fig. HI The data in Fig. H] were obtained for mean 
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Figure 4: (Colour on line) Invariant distribution p versus 9. The legend 
shows the value of x corresponding to each line. 



field U = 8 and disorder strengths y(uf) = a/ ( A^) = 0.02. 

The assumption of a uniform or slightly modulated distribution must 
therefore be radically dropped and this entails that one can no longer neglect 
the logarithm of the ratio (fl4"|) in the expression (fT3|) of the inverse localisation 
length. Taking into account Eq. ( )26|) . in the neighbourhood of the band edge 
the inverse localisation length (Tl3|) can be reduced to the form 



A ~ -( log 

a 



D r 



sm(9 n+1/ 



sin 



(33) 



On the other hand, using the angular map in Eq. (1271) . one can write the sine 
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ratio in Eq. ( 13"3"j) as 



sm {V n+ i, 
sin (0 n ) 



1 + e cot (8 n ) H — - sin (6* n ) cos (6 n ) 



a 



sm 



cos 



sm 

2 



(34) 



+ • 



Substituting the approximate identities ( 1281) and ( 1341) in Eq. ( 1331) . one obtains 



A ~ -( cot 
a 



The average can now be computed using the invariant distribution 
final result is [12] 



the 



£ J 


POO 

1 x 1 / 2 exp 




-2x ( 




dx 


2a i 
J 


poo 

x~ x l 2 exp 


-2x 




dx 



(35) 



Away from the band-edge on the length scale set by the disorder strength, 
i.e., for k — > oo, Eq. (135]) reduces to 



A 



(7 



8ae 2 



which coincides with the form of the standard expression f|T8l) in the limit 
ka = e — > + . On the other hand, close to the band-edge (on the length scale 
defined by the disorder strength), i.e., for x — > 0, Eq. ( 1351) gives 



A 



6 V3 



2ar (1/6) 



2\l/3 



which exhibits the same anomalous scaling found in the Anderson model at 
the band-edge [15], [12]. This correspondence is a consequence of the fact that 
in both models at the band edge the invariant distribution for the angular 
variable has the form ( 132|) and the ratio ip n+ i/ip n reduces to the same function 
of 6. 
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4 Existence of anomalously localised states in 
the Kronig-Penney model 

In this section we discuss how specific cross-correlations between the two 
kinds of disorder (structural and compositional) can endow the Kronig- 
Penney model ([T]) with electronic states whose amplitude, away from the 
localisation centre no, decays like a stretched exponential. More precisely, 

one has 

\ifi n \ ~ exp (-D^\n-n \) (36) 

where D is a constant. This corresponds to a stretched exponential exp(— \x\ a ) 
with stretching exponent a = 1/2. The phenomenon has its counterpart in 
the band-centre anomaly which occurs in the Anderson model with purely 
off-diagonal disorder [15] . 

As remarked in the previous section, the Kronig-Penney model ([1]) has 
an equivalent tight-binding model. The correspondence is easily established 
by eliminating the momenta from the map ([5]); in this way, with the obvious 
substitution x n — > ip n , one obtains the equation 



rVWl + r~~7~~i~~A vrVVi-l 



sin [q (a + A n )] sin [q (a + A re _i)] 

cot [q (a + A n )] + cot [q (a + A n _i)] + - ^ — \ 



(37) 



It is convenient to express the coefficients of Eq. ( 1371) as sums of their mean 
values and of fluctuating terms with zero average. Eq. ( 1371) then assumes the 
form 

(1 + J n ) VWl + (1 + 7„-l) VVi-l + ^n^n = Ei[) n (38) 

with being a deterministic function of the wavevector q defined by the 
identity 

U/q + 2( cot [g (a + A n )]\ 

E(g) = — ^ ^ (39) 

1/sin [q (a + A n )] 



while the symbols 7 n and e n stand for the energy- dependent random variables 

1/sin [q (a + A n )] 
1/ sin [q (a + A n )] 



7n(g) = iZ^lg>±M_ _ ! (40) 
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and 



s n (q) = ~, r {2(cot[g(a + A„)]\-cot[g(a + A n )] 

/l/sin^Ca + A^] X 1 x " 1 



- cot [q (a + A n _i)] - — 

(41) 

Eq. f|55|) shows that the tight-binding counterpart of the Kronig-Penney 
model ([1]) is an Anderson model with both diagonal and off-diagonal disorder. 
Note that, in the absence of structural disorder, the random variables ([40]) 
vanish and the Kronig-Penney' s analogue becomes the ordinary Anderson 
model with diagonal disorder. For purely compositional disorder, therefore, 
one could have predicted a priori the existence of the anomalies discussed in 
Sec. El 

We now focus our attention on the case in which the compositional dis- 
order has the form 



u n = 2q c (^ cot [q c (a + A n )] J - q c cot [q c (a + A n )] - q c cot [q c (a + A n _i)] 

(42) 

where q c represents the wavevector fulfilling the condition 

E(q c ) = 0. (43) 

We stress that condition (T42"]) introduces special cross- correlations between 
the two kinds of disorder in the Kronig-Penney model (pQ). We also remark 
that the weak-disorder condition (T5]) ensures that a solution of Eq. fT45]) exists. 
In fact, taking into account the band-structure relation ([8]), within the limits 
of the second-order approximation one can cast Eq. (143]) in the form 

q c (Al)U 
cos (k c a) ~ — . 

2 sin(g c a) 

If the structural disorder is weak enough, the right-hand side of this equation 
is less than one and this implies that a Bloch vector k c exists such that q(k c ) 
is the solution of Eq. (1431) . Actually, a perturbative calculation shows that 
the Bloch vector k c in the positive half of the first Brillouin zone is 

^•-3!®!> (44) 
2a 2a sm(ga) 
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with q = q (f-) . For weak disorder the deviation of k c from 7r/2a is not large 
and therefore q 2 — q 2 , which is close to centre of the energy band. 

If the compositional disorder has the form (j42p . it is easy to see that, when 
the electron energy takes the critical value q 2 identified by the condition fj43|) . 
Eq. fl38|) becomes 

[1 + 7fi(?c)] VWl + [1 + 7fi-l(ft)] V>n-1 = 0, (45) 

which has the same form of the Schrodinger equation for the Anderson model 
with purely off-diagonal disorder and zero energy. For zero energy, the lat- 
ter model is known to have an electronic state which exhibits anomalous 
localisation, because it is localised but decays away from the localisation 
centre n according to Eq. (pjj) [T9] . We stress that the identity of Eq. ( 14"5|) 
with the zero-energy Schrodinger equation for the Anderson model with only 
off-diagonal disorder ensures that, when the compositional disorder has the 
form (j42p and the energy takes the critical value q 2 , the Kronig-Penney model 
also has an anomalously localised state whose amplitude decays exponentially 
with the square root of the distance from the localisation centre. 

This property can be heuristically justified with the observation that 
Eq. (l4"5j) implies that 

n-l 

log \ip2n\ = log \ipo\ + ^2 t lo S I 1 + 72«(?c)| - log |1 + 7a+i(g c )|] (46) 

1=0 

By invoking the central limit theorem, one can therefore conclude that, for 
large values of n, the random variable log \ip n \ has zero average and a variance 
which increases linearly with n. We stress that the previous argument holds 
even if the disorder exhibits long-range self-correlations of the form 

1 1 

x(l) = ; [ sin (ttc 2 - sin (ttc-J)] (47) 

C 2 — C\ TIL 

where C\ and c 2 are real numbers such that < C\ < c 2 < 1. In fact, 
weaker forms of the central limit theorem can be applied to sums of correlated 
random variables, provided that the correlations decay fast enough [24, 25J. 
Specifically, given a succession of zero-average, correlated random variables 
{x n }, let Sn = J2n=i x n be the sum of the first iV terms of the succession. 
The minimal condition for the mean square of Sn to grow linearly with 
N, i.e., ((Sn) 2 ) ~ N, is that the power spectrum of the succession {x n } 
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be finite at the origin. The power spectrum corresponding to the binary 
correlator (1471) is 



W{ka) 



c 2 - Ci 




if ka G 



" 7T 


7T" 


U 


7T 


7T" 


Ci — 


C2 2. 


7T — Co — 


7T — Ci — 


L 2 






L 2 


2J 



otherwise 



(48) 



and vanishes at the origin; one can therefore conclude from Eq. (j4U 

log | ^„ | - 



that 



n 



even if the disorder is correlated. 

The numerical data confirm the conclusion that the Kronig- Penney model (pQ) 
has an anomalously localised state for q = q c when the compositional dis- 
order takes the special form (14"2l) . This can be seen from Figs. [5] and El which 
show how the Lyapunov exponent vanishes when the Bloch wavevector takes 
the critical value k c . The difference between Figs. and lies in the fact 
that in the case corresponding to Fig. [5] the structural disorder is not self- 
correlated, while the data represented in Fig. [6] were obtained for structural 
disorder with long-range self-correlations of the form (T4T1) with c\ = 3/10 
and c 2 = 1. 

All the numerical data presented in this section were obtained for a 
Kronig-Penney model with mean field U = 4 and structural disorder charac- 
terised by a uniform distribution 



l/W if A n g [-W/2, W/2] 
otherwise 



(49) 



with width W = 0.1732 . . ., corresponding to a disorder strength a/ (A^) = 
0.05. For the above-specified values of the mean field and of the disorder 
strength, formula (jUJ) gives a value of the critical Bloch wavevector ap- 
proximately equal to k c ~ 0.49527r/a, in relatively good agreement with the 
numerically obtained value k c ~ 0.49967r/a. 

When the box distribution ( 149]) is chosen for the displacements A n , the 
function (1391) becomes 



UW + 2 log 



E(q) 



sin [q (a + W/2)} 
sin [q (a - W/2)} 



sin(qa) + sin(qW/2) 



(50) 



log 



sin(ga) — sm(qW/2) 
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71/10 n/5 3tt/10 2tz/5 n/2 3n/5 InllO An/5 9n/10 

ka 

Figure 5: (Colour on line) Inverse localisation length A versus ka for the 
case of structural disorder without self-correlations. The dashed line (1) 
represents numerical data, while the solid line (2) corresponds to Eq. fl54|) . 
The inset shows the anomaly in greater detail. 



while the random variables ( |40l) and ( I4TI) take the forms 

, s i qw 

ln{q> sin [q (a + A n )l , sm(qa) + sm(qW/2) { > 

log ■ 



sin(ga) — sin(gH / /2) 



and 



e iq) = I L sin[g(a + ^/2)] 

n[q) sin(ga) + sm(qW/2) \ ° 8 sin [q (a - W/2)] 

° g sin(ga) - sin(qW/2) (52) 

- qW cot [q (a + A n )] - qW cot [q (a + Vi)] ~ Wu r 



We would like to stress that, when performing numerical calculations, one 
should work with the exact form fl50|) of E(q) and the exact expressions fj5T|) 
and fl52|) of the coefficients 7„ and e n , even if the disorder is weak. Obvi- 
ously, the explicit expressions of these magnitudes depend on the distribution 
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7Z/5 371/10 27T/5 7Z/2 37T/5 77T/10 47T/5 

Figure 6: (Colour on line) Inverse localisation length A versus ka for the 
case of self-correlated structural disorder. The dashed line (1) represents 
numerical data, while the solid line (2) corresponds to Eq. f )54|) . The inset 
shows the anomaly in greater detail. 



chosen for the variables A n and must be modified if the box distribution is 
replaced with another one. Whatever distribution is adopted, however, it 
is important that the corresponding exact expressions of E(q), 7„ and e n 
be used. Using second-order approximations for E(q) and the random coef- 
ficients fj4"T|) works relatively well for most values of the energy, but fails 
at the critical point, because the neglected higher-order corrections produce 
non-zero diagonal terms in Eq. (1381) which, in spite of being very small, pre- 
vent the electronic state from being anomalously localised and the Lyapunov 
exponent from vanishing completely. 

As a side remark, we would like to add that the conditions (TJ2]) and fj43l 
ensure the existence of an anomalously localised state also for disorder of ar- 
bitrary strength, provided that Eq. ( ]4"3l) has a solution q c inside the allowed 
energy bands. This is confirmed by numerical calculations (which, incident- 
ally, also show that when disorder is not weak the critical value of the energy 
need not be close to the band centre). To determine the conditions which 
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guarantee the existence of such a critical value of the energy in the general 
case is not an easy task, however; for the sake of simplicity, we therefore re- 
strict our attention to the case of weak disorder, for which the band structure 
is approximately given by Eq. (JSJ) and Eq. ( |4"3]) does have a solution. 

In both Fig. [5] and [6] the numerical data are compared with the theoretical 
predictions derived from the general result ffTTl) . Because in this formula 
terms of order higher than the second are neglected, in its evaluation we 
replaced the exact expression (1421) of the compositional disorder with its 
second-order approximation 

Q 2 

u n ~ , 2 7 - (A„ + An_i) . (53) 
sin (q c a) 

In passing, we observe that Eq. (|53|) implies that the cross-correlator X3(0 is 
not an even function of its argument. When the compositional disorder has 
the form (1531) . the inverse localisation length (IT7|) becomes 



8a 



l-2/^cos(M 
Uqsm {q c a) 



W 2 (ka). (54) 



The power spectrum ^(fca) in the previous formula reduces to W2(ka) = 1 
in the case of uncorrelated structural disorder, while it is of the form (I4"8l) 
in the case of self-correlated structural disorder represented in Fig. [6j Note 
that in the latter case the long-range correlations (1471) create two mobility 
edges at ka = 37r/10 and ka = 7n/10. As can be seen from Figs. [5] and EJ 
the theoretical formula ( 1541) works reasonably well everywhere, except in a 
small neighbourhood of the critical value k c . This failure must be ascribed 
to the fact that Eq. (1541) . as its parent expression (fTTj) . is not valid for Bloch 
wavevectors lying close to the rational value 7r/2a, where resonance effects 
play a non-negligible role. 

As a last comment on the specific features of the localisation length in 
the special case in which the compositional disorder is related to the struc- 
tural one by Eq. (|42|) . we observe that Eq. (I54T) predicts the existence of a 
delocalised state for the Bloch wavevector k* identified by the condition 

2 ^in( 9(fe >) o 
q{k*) sin {q c a) 

This is confirmed by the numerical computations, as can be seen from both 
Figs. [5] and El 
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Figure 7: (Colour on line) (log 1^1) versus n. The dashed line corresponds 
to the critical value of the Bloch wavevector, k = k c ~ 0.49967T, while the 
dotted line corresponds to a Bloch wavevector k = 0.50367T. 



To conclude this section, we would like to add some numerical evidence of 
the stretched-exponential behaviour of the tails of the anomalously localised 
state. We have numerically solved Eqs. ( 138|) and ( j4"5i) as initial-value prob- 
lems; this corresponds to constructing the electronic states with the transfer- 
matrix technique. In the anomalous case one expects log \ip n \ to behave as 
the position of a random walker, 1.6. , cLS cL random variable with constant zero 
average and a second moment linearly increasing with n. For k ^ k c , on the 
other hand, the solution of Eq. (138]) should behave as \ip n \ ~ exp(An), leading 
to an increase with n l of the Z-th moment of log \ ip n \. This is confirmed by the 
numerical data for the first two moments of the variable log \ ip n \, represented 
in Figs. [7] and |HJ In both Figs. [7] and [HJ we considered the behaviour of 
log \ip n \ as a function of n for two Bloch wavevectors, i.e., the critical vector 
k = k c ~ 0.49967T and the vector k = 0.50367T, which is close to the critical 
value but not identical to it. The moments of log \ip n \ were computed with 
an average over 1000 disorder realisations. For the sake of simplicity, we con- 
sidered compositional disorder without self-correlations. Both the first and 
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Figure 8: (Colour on line) (log 2 \ip n \) versus n. The dashed line corresponds 
to the critical value of the Bloch wavevector, k = k c ~ 0.49967T, while the 
dotted line in the inset corresponds to a Bloch wavevector k = 0.50367T. 



the second moment of log \ip n \ behave as expected, corroborating the conclu- 
sion that the tails of the electronic state at the critical point are described 
by Eq. (J36]). 



5 Conclusions 

In this work we analyse the anomalous behaviour of specific electronic states 
in the Kronig-Penney model with weak compositional and structural dis- 
order. In every case we discuss the analogies with the corresponding phe- 
nomena in the Anderson model. 

We first show that the localisation length deviates from the prediction of 
the standard formula (TIT]) when the Bloch vector assumes the value k = n/2a, 
which corresponds to an energy close to the band centre. This discrepancy 
is due to the same resonance effect which occurs in the standard Anderson 
model at the band centre; in both models this effect produces a modulation of 
period 7r/2 of the invariant distribution of the angle variable of the associated 
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Hamiltonian map. This modulation leads to the band-centre anomaly of the 
localisation length. 

We make use of the Hamiltonian map approach also to analyse the elec- 
tronic states at the band edge. We find again that the same anomalous 
behaviour originally found in the Anderson model is present in the Kronig- 
Penney model. In both systems the most relevant feature at the band edge is 
the anomalous scaling of the localisation length with the disorder strength. 

We finally use the correspondence between the Kronig- Penney model and 
the Anderson model with diagonal and off-diagonal disorder to conclude that 
in the former system specific cross-correlations between the two kinds of dis- 
order generate an electronic state, close to the band centre, whose tails decay 
as stretched exponentials. This state is the analogue of the anomalously loc- 
alised state which occurs at the band-centre in the Anderson model with 
purely off-diagonal disorder. 

J. C. H.-H. and L.T. gratefully acknowledge the support of the CONACyT 
grant No. 84604 and of the CIC-2009 grant (Universidad Michoacana). The 
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